Use the interactive map to check population changes.
points_bau4326 <- points_bau |>
#filter(BASIN_N == "Orinoco") |>
arrange(BASIN_N, SUBBASI, BB_ID, REACH_ID) |>
filter(row_number() %% 5 == 1) |> st_as_sf() |> st_transform(4326)
points_bau4326_low <- points_bau4326
# colour palette
pal <- colorFactor(
palette = c("#7274C1", "#A3720E"),
domain = points_bau4326$flag_50_42y
)
# interactive map
leaflet::leaflet(data = points_bau4326) |>
addTiles() |>
addCircleMarkers(color = ~pal(flag_50_42y),
stroke = FALSE, fillOpacity = 0.4,
clusterOptions = markerClusterOptions(),
group = "points_bau4326"
) |>
groupOptions("points_bau4326", zoomLevels = 1:8) |>
addCircles(data = points_bau4326_low,
color = ~pal(flag_50_42y),
group = "points_bau4326_low"
) |>
groupOptions("points_bau4326_low", zoomLevels = 9:15)